Method to Look for Imprints of Ultrahigh Energy Nuclei Sources 
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We propose a new method to search for heavy nuclei sources, on top of background, in the Ultra- 
High Energy Cosmic Ray data. We apply this method to the 69 events recently published by the 
Pierre Auger Collaboration [l[ and find a tail of events for which it reconstructs the source at a few 
degrees from the Virgo galaxy cluster. The reconstructed source is located at ~ 8.5° from M87. 
The probability to have such a cluster of events in some random background and reconstruct the 
source position in any direction of the sky is about 7 x 1CP 3 . The probability to reconstruct the 
source at less than 10° from M87 in a data set already containing such a cluster of events is about 
4 x 1(T 3 . This may be a hint at the Virgo cluster as a bright ultra-high energy nuclei source. We 
investigate the ability of current and future experiments to validate or rule out this possibility, and 
discuss several alternative solutions which could explain the existing anisotropy in the Auger data. 
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I. INTRODUCTION 

The HiRcs and Auger data have proved the existence of 
a cutoff at the highest energies in the Ultra-High Energy 
Cosmic Ray (UHECR) spectrum @, [|. This indicates 
that UHECR sources are astrophysical. Due to the GZK 
cutoff or photo-disintegration, they must be located in 
the local Universe (at distances r < 100 Mpc). These 
astrophysical sources, which are still unknown, should 
be located in the Large Scale Structure (LSS) of matter. 

The Pierre Auger Collaboration recently reported a 
shift towards a heavier composition in the Ultra- High En- 
ergy Cosmic Ray spectrum at the highest energies, above 
a few times 10 19 eV Q. The analysis of the Yakutsk EAS 
Array muon data is also in agreement with this observa- 
tion However, these results are still controversial: 
the measurements of HiRcs experiment ml and some pre- 
liminary studies of the Telescope Array [7| are consistent 
with a proton composition. 

For the moment, methods to search for the sources of 
UHECRs have been presented forproton or light nuclei 
primaries. Sec for example Refs. 

In case UHECR are heavy nuclei, looking for their 
sources would be a harder task: for example, 60EeV 
iron nuclei behave as ~ 2EeV protons in the Galac- 
tic Magnetic Field (GMF), due to their similar rigidities 
E/Z. Refs. [HUH] studied the propagation in the GMF 
of particles with such low rigidities, while Refs. [HI, HH 
discussed the effect of varying the UHECR composition 
on the correlation of Auger events with active galactic 
nuclei. 

The GMF displays both a large scale and a random 
small scale structure, which are respectively known as 
the regular and the turbulent components. 
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In Ref. [l7[, T. Stanev suggested one of the first mod- 
els of the regular GMF, describing analytically with 
logarithmic spirals the field structure in the Galac- 
tic disk. Other models of the disk field were formu- 
lated by D. Harari et al. and by P. Tinyakov and 
I. Tkachev [3. Then, M. Prouza and R. Smida built a 
model which adds a halo contribution made of toroidal 
and poloidal fields [19|,|20|. The implications on the GMF 
modeling of recent rotation measure maps were reported 
by J. L. Han et al. [21H231 ] . Other models have been sug- 
gested for the disk field, such as a toroidal field consisting 
of concentric rings by J. P. Vallee et al. [2411 and another 
axisymmctric field by L. Page et al. j25j . Some recent 
models for the disk field also display a spiral pattern 
based on the structure of the NE2001 thermal electron 
density model (J. C. Brown et al. J26[) or on the spiral 
structure of the Milky Way (Y. Y. Jiang et al. H3) de- 
duced from HII regions and giant molecular clouds (28| . 
X. H. Sun et al. proposed several GMF models and con- 
fronted them with the data [2t|[3(|. However, currently 
no theoretical GMF model can reasonably well fit all ex- 
perimental data, as Refs. [U [HJ show. 

References [33l . |34| discuss the turbulent component 
modeling and Ref. [35| investigates its spectrum. Its im- 
plications on the propagation of UHECRs are studied in 

Refs. jHSIalls^- 

In many cases, one may not detect the signatures from 
Ultra-High Energy (UHE) nuclei sources without a more 
precise knowledge on the Galactic Magnetic Field than 
currently available (l4j | . Astronomy with UHE heavy nu- 
clei can look very different from astronomy with light pri- 
maries. In particular, multiple images of the same source 
can appear even at the highest energies. The images 
of nearby galaxy clusters in some recent GMF models 
have been shown in Ref. [lH, illustrating the challenges 
of heavy nuclei astronomy. 

However, the images of some UHE nuclei sources could 
be detectable in favorable cases, without requiring an im- 
proved knowledge of the GMF. We propose in this paper 
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a method to look for heavy nuclei sources, in such situ- 
ations. It is an extension to the case of heavy nuclei of 
the method we presented in Ref. @ for protons and light 
nuclei. 

We apply the method to the list of 69 events with ener- 
gies E > 55 EeV recently published by the Pierre Auger 
Collaboration [l| . We find in this data a cluster of events 
for which the reconstructed source lies near Virgo, which 
is in line with the supposition of Ref. (37} ■ We generate 
sets of 69 background events following the exposure and 
spectrum of Auger. The probability to have such a clus- 
ter of events and reconstruct the source position in any 
direction of the sky is about 7 x 1CP 3 . Assuming that 
the cluster already exists in the data due to another rea- 
son, we study how often the method would reconstruct 
the source at less than 10° from M87. This probability 
is about 4 x 10~ 3 . The combined probability of having 
such a cluster and reconstructing the source position at 
less than 10° from M87 is about 3 x 10 -5 . Being the 
largest galaxy cluster in the local universe and hosting 
the powerful active galaxy M87, the Virgo cluster is, the- 
oretically, a good candidate for being home of one or 
several source(s) of UHECR. 

Nevertheless, both the "limited" amount of data and 
the poor knowledge of the GMF global geometry prevent 
us to conclude firmly whether the detected tail of events 
really comes from Virgo or not. We analyse the ability 
of current and next generation experiments to test this 
possibility. We also review alternative solutions which 
could explain the Auger data. 

This paper is structured as follows: In Section UIl we 
present a method to search for detectable imprints of 
UHE nuclei sources. We apply it to the recently pub- 
lished Auger data, in Section Mil Wc find with this 
method a signal which may hint at the Virgo cluster as 
a bright UHE nuclei source. This result is tested with a 
"blind-like" analysis in Section [TV] by analysing succes- 
sively the lists of 27 and 69 — 27 = 42 events recorded 
by Auger. In Section |Vj we present a detailed and criti- 
cal analysis of the possibility that Virgo may be a UHE 
nuclei source. We also discuss alternative reasons which 
could explain the anisotropy in the Auger data. 



II. METHOD TO SEARCH FOR HEAVY 
NUCLEI SOURCES 

At the highest energies, for a source of UHE protons or 
light nuclei located far enough from the Galactic plane, 
the GMF approximately shifts its events in a sector- 
shaped region on the celestial sphere, on one side of the 
source The vertex of the sector is, theoretically, lo- 
cated at the source position, and its opening angle de- 
pends on the ratio of deflections due to the turbulent 
and regular components of the GMF. In a first approxi- 
mation, the angular distance between the source and an 
event of energy E is proportional to 1/E Q . The associ- 
ated proportionality factor will be called the "deflection 
power" of the regular GMF, V, in the following. 
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FIG. 1. Arrival directions of ultra-high energy iron nuclei 
emitted by the Hydra galaxy cluster, in Galactic coordinates, 
and deflected in the PS regular GMF model. Colors corre- 
spond to the energies of the emitted nuclei: dark blue stands 
for 60 EeV nuclei, light blue: 70 EeV, green: 80 EeV, yel- 
low: 90 EeV, orange: 100 EeV, red: 120 EeV, and magenta: 
140 EeV. The Hydra cluster is represented by the black disk. 
The image which is approximately "enlarged proton-like" is 
surrounded by the red box. 

Contrary to the naive idea that sources of heavy nu- 
clei would display the same features only enlarged by a 
factor Z, Refs. [10l - [l2l . (3] point out that their images 
often have more complicated shapes, even at energies 
E > 60EcV. For instance, sources can have several im- 
ages. They can appear above or below a given energy 
threshold, and merge into one single image when the en- 
ergy is increased. They can also be strongly distorted 
on the celestial sphere and display an energy ordering 
far from the 1 /E behavior expected close to the ballistic 
regime -see details in Ref. Moreover, these images 

are very dependant on the considered model of GMF. 
Hence, for nuclei sources, a better knowledge than cur- 
rently available of the GMF geometry would be needed 
in many cases, in order to find out an efficient and par- 
ticular algorithm to detect their events and reconstruct 
their positions. 

Meanwhile, one can still try to find simple proton- 
like orderings, with ~ Z times larger angular scales. 
As shown below, this can indeed happen in some fa- 
vorable specific cases. For some types of GMF models 
and for some positions on the sky, at least one image of 
an iron source can look like a more or less roughly en- 
larged proton-like image at high enough energies (E > 
50 — 60 EeV). We checked this point by computing the 
iron images of nearby galaxy clusters in the three re- 
cent models of the regular GMF which are considered in 
Ref. [lH ■ One example of roughly "enlarged proton-like" 
image is the main ima ge o f Hydra cluster in the Prouza 
and Smida (PS) model [19l |2C| . shown in Fig.[TJ This im- 
age is surrounded by a red box. The hydra galaxy cluster 
is represented by the black disk, and colors correspond to 
the energies of the emitted nuclei: dark blue stands for 
60 EeV nuclei, light blue: 70 EeV, green: 80 EeV, yellow: 
90 EeV, orange: 100 EeV, red: 120 EeV, and magenta: 
140 EeV. 
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For this paper, we will focus on the favorable case of 
approximately "enlarged proton-like" images of sources. 
We leave the more frequent but more complicated cases 
for future works. 

The method we proposed to look for proton and light 
nuclei sources is presented in details in Section 3.1 of 
Ref. @. In this work, we slightly modify this method in 
order to optimize it to the search for heavy nuclei sources. 
Since we want to scan over all the free parameters of the 
method, we try to have as few parameters as reasonably 
possible. There are four of them. A schematic image of 
this method is drawn in Fig. [21 The source S is repre- 
sented by a black disk. 

Our procedure starts selecting an event with energy 
Ei > 10 20 eV. We will call it the "highest energy event" . 
It is denoted by "1" in Fig. H 

Let us consider such an event. We do an assumption 
on the typical value of the local regular GMF deflection 
power, T>, and only consider the events which angular 
distance to the highest energy event, d, satisfies 



d< R = 



V 



55EcV 



V 
~E~i 



(1) 



The next step is to search for the events which energy E2 
and distance d also fulfill the condition 



, V V 

d< 

E 2 Ei 



(2) 



The events satisfying this latter condition are tested one 
after another, by decreasing energy order, with the pro- 
cedure described below. 

Let us start with the event which has the highest en- 
ergy among them. It is denoted by "2" in Fig. [2] 

In the following, we focus on the events located in a 
given sector-shaped region of the sky. Its direction is 
defined by this event "2". This region is highlighted in 
grey in Fig. [2] We define it as an extension to spherical 
geometry, of a circular sector which vertex is located at 
the position of the highest energy event ("1" in Fig. [2]). 
More precisely, such region is the sub-region of a spher- 
ical lune with the highest energy event on one vertex. 
This sub-region contains the points of the spherical lune 
which angular distance to this vertex satisfy Eqn. (flj. 
We will refer to this region of the sky as the "sector" in 
this paper. Its opening angle is given by the second free 
parameter of the method, 9, and its extension by the 
angular distance R -see Eqn. (fTJ). Its central axis, which 
divides the opening angle in two equal parts, is defined 
by the line containing both the events "1" and "2" . 

Let us define the correlation coefficient Corr(X' , 1/E) 
for the events in the sector, where E denotes their ener- 
gies and X' their angular distances to the vertex of the 
sector. X' is represented by the red axis in Fig. [21 The 
two last free parameters of the method are the minimal 
number M of events in the sector and the minimal value 
Cmin of the correlation coefficient: If there are more than 
J\f events in the sector, and that CorriX' , 1/E) > C m i n , 
there is a detection. Otherwise, the second event in the 
ordered list of events satisfying Eqn. ([2]) is tested. The 




FIG. 2. Example of a UHE heavy nuclei source detection, in 
case of an "enlarged proton-like" image. The extended source 
S is shown by the black disk, and its cosmic rays are repre- 
sented by decreasing energy order, by the magenta filled and 
red open boxes, green open circles and blue crosses. The "sec- 
tor" , with an opening angle and extension R, is highlighted 
in grey -see text for details. "1" and "2" respectively point at 
the highest energy event and at the selected event in the list 
satisfying Eqn. ([2]). X' axis is drawn in red. S' denotes the 
reconstructed source. 



procedure continues until either there is one detection, 
or all the events in the list are tested. In the latter case, 
there is no detection. 

In case of detection, we reconstruct the source position 
as depicted in Section 3.3 of Ref. @. The source is re- 
constructed along the X' axis. This axis is not the exact 
central line of the sector. It is the axis which contains 
both the center of mass of all cosmic rays in the sector 
and the vertex of the sector. The position of the recon- 
structed source is given by the fit of 1/E versus X' . It is 
represented by a thick red cross, S', in Fig. [21 

For heavy nuclei with energies E > 55 EcV, deflections 
on the celestial sphere can easily reach several tens of 
degrees. Even in good cases, this often leads to strong 
deviations to the linear shape of images and to the 1 /E 
ordering of events. That is why, one should not expect 
for heavy nuclei sources the same excellent precision on 
the reconstruction of the source position as for proton 
sources. 

To summarize, the method used in this paper has 4 
free parameters: 

1. The deflection power T>, and the opening angle 
O. The best value for D is mostly related to the 
strength of the regular GMF. The best value for 
O mostly depends on the ratio of deflections in 
the turbulent and regular components of the GMF. 
While the values of these contributions are not pre- 
cisely known due to the lack of knowledge on the 
GMF, their most probable ranges can be inferred 
from the literature. 

2. The minimum number of events in the sector, A/", 
and the minimum value of the correlation coeffi- 
cient Corr(l/E, X') for the events in the sector, 
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Cmin- Below these values, the considered features 
are rejected by the method. 



III. APPLICATION TO THE DATA OF THE 
PIERRE AUGER OBSERVATORY 

The Pierre Auger Collaboration recently released in 
Ref. [l| a list of 69 events with energies higher than 
55 EcV. They were recorded with a total integrated ex- 
posure of 20,370 km 2 -sr-yr. In Fig. [3] (upper panel), we 
plot in Galactic coordinates the positions of the first 27 
events. They correspond to the first released data set of 
Refs. [38l. [39i|. renormalized as in Q. We plot in Fig. 
(lower panel), all the 69 events (current data set). Events 
with energies E > 10 20 eV, 10 199 eV< E < 10 20 eV, 
10 198 eV< E < 10 19 ' 9 eV and 55 EeV < E < 10 19 8 eV 
are respectively represented by filled magenta boxes, red 
open boxes, green open circles and blue crosses. 

The most visible feature in the Auger data is an over- 
density of events in the region —60° < I < —30° and 
0° is b ^ 30°. It was first discussed for the data set of 
27 events in Ref. [4(| ■ For the 69 events data set, it was 
studied in Ref. In the following, we will call this part 
of the sky the "Cen A region" . An important point was 
noted in Ref. [37j ■ It shows that if one excludes this over- 
density, the rest of the sky could still be compatible with 
isotropy. The significance of the overdensity was com- 
puted both in Refs. and H3- Ref. [1 found 4% with 
a Kolmogorov-Smirnov test, and Ref. [331 3% (resp. 2%) 
with 3-point (resp. 4-point) autocorrelation functions. 

We shall now analyse the set of 69 events with the 
method presented in the previous section. We scan the 
Auger data over all possible combinations of the 4 pa- 
rameters of the model. We use discretized sets of values 
for each parameter. 

We take N G {4, 5, 69}. For the minimum value 
of the correlation coefficient, we take a set with steps 
equally spaced by 0.1: C min G {-1, -0.9, 0.9}. Ex- 
isting measurements of the GMF allow to give typi- 
cal estimates of ~ 1 — 2.5° deflections on the sky for 
10 20 eV protons [2(|. Then, assuming the heavy nuclei 
primaries to be iron (Z = 26), we take for this analy- 
sis V G {26°, 39° 52°, 65°} x 10 20 eV. According to the 



results of Ref. [34| on the relative contributions of the 



regular and turbulent components to the UHECR de- 
flections, < 80° should be sufficient. Then, we take 
6 G {10°, 20°, ...,80°}. 

We confront below the data with 4.7 x 10 7 Monte Carlo 
simulations of random background. This background 
is made of exactly 69 events and follows the exposure 
of Auger experiment, i.e. its local density statistically 
follows the exposure. Its energy spectrum follows the 
spectrum of the 69 Auger events, distributed in four 
logarithmically spaced bins: 55 EeV-10 198 eV, 10 19 8 eV- 
10 199 eV, 10 199 eV-10 20 eV and above 10 20 eV. The ex- 
act distribution of energies within each of the first three 
bins does not significantly change the results below. The 
spectrum above 10 20 eV is poorly known, and Auger has 
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FIG. 3. Images, in Galactic coordinates, of Auger events with 
energies E > 55EeV, published in Ref. [l[. Upper panel: 
First data set of 27 events [H, [3S|, with the renormalized 
coordinates of Ref. Lower panel: Full data set of 69 
events. Black disk for the position on the sky of the Virgo 
cluster. Filled magenta boxes stand for events with energies 
E > 10 20 eV, red open boxes for 10 199 eV < E < 10 20 eV, 



green open circles for 10 19 
crosses for 55 EeV < E < 



eV< E < 10 1 



eV and blue 



10 



10.8 



eV. The red arrows show 



the direction along which events emitted by Virgo would be 
shifted in the GMF if Cen A region events are nuclei from 
Virgo. 



only recorded 3 events at such energies. We take here 
for the events in the bin E > 10 20 eV a E~ 4 3 spectrum. 
This spectrum was proposed in Ref. 4l| for events with 
E > 10 19 5 eV. The value for the maximum observable 
energy, E max , is chosen according to results on propaga- 
tion of nuclei. The figure 1 of Ref. [42| shows that the 



iron nuclei propagation length rapidly falls below a few 
Mpc for energies above 3 x 10 20 eV. Therefore, we take 
E max = 10 20 - 5 ~ 3 x 10 20 cV. We checked that taking 
lower values for E ma , x would only increase the significance 
of the signal detected below. 

For one of the three events with energies above 10 20 eV, 
we find an interesting signal which is shown in Fig. 01 The 
coordinates of this event are : E = 142 EeV, I = —57. ,2°, 
b = 41.8°. It is located at ~ 34° from the center of 
Virgo. It plays the role of the "highest energy event" 
in the method. The best configuration is obtained for 
N = 13, C min = 0.6, V = 39° x 10 20 eV and 9 = 40°. 
Among all possible configurations in the data, it is the 
one which has the lowest probability to be reproduced by 
the background. The value for Corr(X' , 1/E) is ~ 0.66. 

We computed the probability to obtain an at least as 
good feature in the background. Out of 4.7 x 10 7 gener- 
ated sky maps of background, we found 3868 of such fea- 
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FIG. 4. Portion of the celestial sphere with the Auger data 
in Galactic coordinates, for the case of the best sector with 
9 = 40° and V = 39° x 10 20 eV. The 13 events contained in 
the sector are surrounded by magenta circles. The blue line 
represents the X' axis and the thick red cross, the position of 
the reconstructed source (I ~ —106°, b ~ 72.5°). The Virgo 
galaxy cluster is denoted by the black disk. Same color code 
for the Auger cosmic rays, as in Fig. [3] 



tures for these 4 fixed parameters. This number has to be 
penalized over all possible values in the ranges of the 4 pa- 
rameters of the method. As pointed out in Ref. [431 ] . one 
cannot know a priori the best values of the scan parame- 
ters. Therefore, one has to take into account any config- 
uration in the background which probability is lower or 
equal to the probability of the specific feature detected 
in the data for N = 13, C min = 0.6, V = 39° x 10 20 eV 
and O = 40°. In this study, for each value of (2?, 9) 
we scan all Monte Carlo realizations of the background 
over all values of (TV, C m - m ). We count all cases for which 
the probability to have a detection with a given value of 
(A/", C m in) is lower or equal to the probability of the best 
sector in the data, 3868/ (4.7 x 10 7 ). After summing over 
all values of (X>, 0), the total number of such cases in the 
background is 311481. With our method, the probability 
of the feature is then: P feature ~ 6.6 x 10~ 3 . 

For this configuration, the reconstructed source is lo- 
cated at ~ 8.5° from M87, at (I ~ -106°, b ~ 72.5°). It 
is near the boundaries of the Virgo cluster, which has an 
apparent radius of ~ 5° on the sky. 

The position of the reconstructed source is drawn from 
the central value of the fit of 1/E versus X' . The uncer- 
tainties due to the fit are of the order of <~ 10°, because 
of the low energy ordering of events in the Cen A re- 
gion. Moreover, even if M87 would be the only source in 
Virgo, the magnetic fields in the cluster can be sufficient 
to significantly deflect trajectories of UHE heavy nuclei 
inside and make shine the whole cluster as an extended 
source [iij . 

The reconstructed position is then compatible with the 
Virgo cluster (or M87) being the source of (most of) the 
considered 13 events. In Fig. SI these 13 events are sur- 
rounded by magenta circles and the reconstructed source 
position is denoted by the thick red cross. 

If we also add the constraint that the reconstructed 



FIG. 5. Renormalized probability P/ < P > to reconstruct 
the source in any direction of the celestial sphere with ran- 
dom background following the Auger exposure and spectrum. 
< P > represents the mean probability P averaged over all 
bins of the sky map. The white circle indicates the 10° region 
around M87. 

source should be located at less than 10° from M87, 
the number of such cases in the background falls to 
15 and 1214, respectively before and after the penal- 
ization. This corresponds to the following probability: 
P/eatMre,d(A/87)<io° - 3 x 1CT 5 . The Virgo galaxy cluster 
was suggested to be a source of UHE nuclei by Ref. [33] • 
The events in the Cen A region and the 142 EeV event 
were supposed to be its image, which is in line with the 
results shown in Fig. 21 

However, one could argue that the overdensity can 
be explained by another reason than events emitted by 
the Virgo cluster. The Cen A region overdensity may 
for example be explained by a magnetic lensing effect, 
and the presence of the nearby event with E > 10 20 eV 
may have triggered artificially the detection. The rel- 
evant value is then the probability to reconstruct the 
source at less than 10° from M87, in case one has already 
such a feature in the data (i.e. at least 13 events, with 
a correlation coefficient > 0.6). This is estimated by: 

Pd(M87)<10° = P feature,d(M87)<10° /P feature — 4 X 10 

Thus, the probability that the source is reconstructed 
near Virgo due to a random fluctuation is ~ 0.4%. 

It is slightly lower than the value it would take in case 
of a random position of the reconstructed source. With 
the background events, the exposure of Auger favours the 
reconstruction of sources within its region of high expo- 
sures. In Fig. [5] we plot the probability P to reconstruct 
the source position in any direction of the sky for random 
background events, renormalized to P/ < P >. < P > is 
the mean probability to reconstruct the source in a given 
bin of the plot, averaged over all directions of the celestial 
sphere. All sectors with a probability of occurrence lower 
or equal to the probability of the best sector in the data 
have been considered. Therefore, the sum of all bins of 
Fig.[5]adds up to J2 P — P feature — 6.6 x 10~ 3 . The part 
of the sky where P > (0.4 — 0.5) x < P > globally corre- 
sponds to the directions in which Auger exposure is non- 
zero. The probability to reconstruct the source in regions 
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FIG. 6. Sector with: Upper panel: 9 — 60°. There are 
15 events in the sector; Lower panel: O = 20°. There are 
8 events in the sector. Same key as in Fig. [4] The source is 
reconstructed at (I ~ -122°, b ~ 74.8°) for 9 = 60°, and at 
(I ~ -78.4°, 6 ~ 72.4°) for 9 = 20°. 



which Auger is blind to is lower, though non-zero. The 
region of maximum probability, P ~ (1.6 — 1.9) x < P >, 
is a circular band within regions of high exposures. The 
white circle surrounds the part of the sky within 10° from 
M87. In most of it, P ~ (0.5 - l)x < P > (blue and 
purple colors). P is slightly larger than < P > in the 
smaller pink subrcgion at lower b and larger I. On av- 
erage, P is slightly lower than < P > in the 10° radius 
region. 

It may be noteworthy to point out the main two other 
possible sector angles 9 that one can consider when 
analysing the Auger data. The value of V = 39° x 10 20 eV 
is left unchanged: 

• = 60° -see Fig. [7| (upper panel): when one con- 
siders this larger sector angle, one selects two more 
events in the overdense region, compared to the 
case 9 = 40°. As visible in Fig. |6l all the cos- 
mic rays which belong to the overdense region are 
taken into account in this configuration. Af = 15 
and C* min = 0.6 (because Corr(X',l/E) ~ 0.61). 
The reconstructed source is located at ~ 11.8° 
from M87, at (I ~ -122°, b ~ 74.8°). The re- 
sults are slightly worse than for 9 = 40° which 
is the "true" minimum for the considered sets of 
parameters, but they are not very far. For this 
case, the probabilities introduced above become: 



FIG. 7. Positions of the reconstructed sources on the celestial 
sphere for the three considered sector opening angles: 9 = 
20, 40, and 60°. Same key as in previous Figures. 



P feature ^ l.lxlO" 2 (respectively 6150 and 503901 
cases out of 4.7 x 10 7 before and after the penaliza- 



tion), P 
2658), P d( M87) 



feature, d(M&7)<10° 



6 x 10" 5 (resp. 30 and 



<io° 



5 x 10~ 3 (0.5%). 



9 = 20° -see Fig. [3| (lower panel): we discuss this 
smaller sector, because the source is reconstructed 
much closer to M87. However, it only takes into ac- 
count 8 points from the overdense region. The re- 
constructed source is located at ~ 2.2° from M87, 
at (I ~ -78.4°, b ~ 72.4°), which is located in 
the Virgo cluster. N = 8 and C m i n = 0.7 (be- 
cause Corr(X',l/E) ~ 0.78). These 8 events be- 
long to a "filamentary" structure which is a denser 
sub-region of the "right" part of the overdense 
region. Knowing if this filamentary structure is 
the real image of the Virgo cluster, instead of the 
whole overdense region, is beyond the scope of what 
one could currently say. We compute the same 
probabilities as above, except that we take into 
account sources reconstructed at distances below 
3° from M87, instead of 10°. Here, P fea ture ^ 
2.1 x 10" 1 (respectively 128469 and 9987251 cases 
out of 4.7 x 10 7 before and after the penalization), 



feature,d(M87)<3° 



9 x 10~ 5 (resp. 53 and 4368), 

Pd(M87)<3° = P feature, d(M87)<3° /P feature — 4 X 

10 -4 . The probability to have in the background 
such a filamentary structure of > 8 events and 
Corr(X',l/E) > 0.7, with V = 39° x 10 20 cV, is 
much higher than the probability to have the fea- 
tures with > 13 or 15 events for 9 = 40° and 60°. 
However, as shows the value of Pd(M87)<3°i once 
one has such a feature, the probability to recon- 
struct the source at less than 3° to M87 is natu- 
rally around 10 times lower than the probability 
to reconstruct it at less than 10° -see the case of 
9 = 40°. 

Fig. [7] shows on the same sky map the three positions 
of the reconstructed sources for the cases: 9 = 20°, 40° 
and 60°. The larger the sector angle, the further from 
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FIG. 8. Generic way to define the second sector which is used 
in the study of Section ITVl Same key as in Fig. [2] R' and 0' 
defined in the text. 



M87 the reconstructed source. This could simply mean 
that the events are not deflected along a straight line on 
the sphere. This would not be surprising for deflections 
of several tens of degrees, even for "enlarged proton-like" 
images [l4j |. On the contrary, this might mean that only 
the events in the sector O = 20° come from Virgo. How- 
ever, having such a thin linear filamentary structure on 
the sky for such large deflections would be hard to realize. 
Let us note that only one event has come in this sector 
in the second Auger data set of 69 — 27 = 42 events. 

No noteworthy feature is found in the data when apply- 
ing the method presented in Section [TT] to the two other 
events with energies above 10 20 eV (115 and 123 EcV 
events). 



IV. CROSS-CHECK WITH A BLIND-LIKE 
ANALYSIS 

We shall now check the result of the previous section 
by doing a "blind-like" analysis. It consists in choos- 
ing the best sector for the first data set released by the 
Pierre Auger Collaboration, which contains 27 events. 
We take the sector for which the measured signal (num- 
ber of events and correlation coefficient Corr(l/E,X')) 
has the lowest probability to be reproduced by the back- 
ground. The method used to define this sector can con- 
tain as many parameters as needed. One fixes the best 
parameters and then does not have to penalize over them. 

The second step consists in analysing the "newer" 69 — 
27 = 42 events of the second Auger data set with this 
fixed "best sector" . 

For proton sources or "enlarged proton-like" images of 
nuclei sources, one can expect that events are roughly 
deflected in a sector which vertex is theoretically located 
at the source position (see Ref. for a full explanation). 
Nevertheless, the source position is a priori unknown. 
Since the highest energy event is near the source position 
for proton sources, one can take the highest energy event 
as the origin of the sector in this case. For UHE heavy 



nuclei deflected in the GMF, the distance between the 
source and its highest energy events is usually estimated 
to be of the order of a few tens of degrees. Taking, as in 
the previous section (see Fig. the highest energy event 
as the origin of the sector gives good results in practice. 
However, one may miss in some cases one or two source 
events at the very border of the sector, notably in the 
regions near the highest energy event. 

In this section, we can take this point into account 
for the method used to define the "best sector" . This 
improved method contains a fifth parameter, named 0' 
below. It starts in the same way as the method of Sec- 
tion IIIII the vertex of a first sector is set on the highest 
energy event and the source position is reconstructed as 
previously. The only difference is that the reconstructed 
source position is now regarded as the vertex of a sec- 
ond sector as shown in Fig. [3] The highest energy event, 
denoted by "1" in Fig. [31 defines the central axis of this 
second sector. It has an opening angle 0' and extends up 
to B! = P/(55EeV). It is highlighted in grey in Fig. |H] 
In this section, it is the sector in which one counts both 
Af and Corr(l/E, X'). Its geometry is more adapted to 
grab the events of an "enlarged proton-like" image of a 
source located near its vertex. 

For the 27 events data set, the best first sector is ob- 
tained for V ~ 41° x 10 20 eV, ~ 40°, and contains 9 
events (including the highest energy event). The corre- 
lation coefficient is Corr(l/E,X') ~ 0.77. The source is 
reconstructed at (I ~ -90.6°, b ~ 71.3°), at ~ 5.3° from 
M87. See upper panel of Fig. |H1 The events in the sector 
are surrounded by magenta circles and the reconstructed 
source position is denoted by the red cross. The prob- 
ability to have with some background at least 9 events 
with Corr(l/E,X') > 0.77, for V = 41° x 10 20 eV and 
= 40°, is Pi - 6 x 10~ 6 . 

The vertex of the second sector coincides with the po- 
sition of the reconstructed source. The best values for 
this sector arc V ~ 41° x 10 20 cV, 0' ~ 30°. For the 
27 events data set, it contains Af = 10 events (including 
the 142 EcV highest energy event) and the correlation 
coefficient is Corr(l/E, X') ~ 0.73. We now fix this sec- 
tor which represents the "best sector" for the first Auger 
data set of 27 events. See middle panel of Fig. [9j where it 
is highlighted in orange. The probability to have at least 
10 events and Corr(l/E, X') > 0.73 in this fixed sector 
with some background made of 27 events, and following 
the Auger exposure and spectrum, is P2 ~ 10~ 8 . 

We can now apply this "best sector" to the second 
Auger data set of 69—27 = 42 events. We find inside Af = 
5 events, with a correlation coefficient Corr(l/E,X') ~ 
0.38. Sec lower panel of Fig. |H1 For the first data set 
of 27 events, the significance was higher than for the 
data set of 42 events. We finally confront this result 
with some random background made of 42 events follow- 
ing the Auger exposure and spectrum. We find that the 
probability to have with the background Af > 5 events 
and Corr(l/E, X') > 0.38 in the "best sector" is equal 
to Pbla ^ 1.8 x 10~ 2 - 2%. 

Let us note that if one adds or removes by hand one 
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FIG. 9. Blind-like analysis conducted with the two consecu- 
tive data sets of Auger (27 and 69 — 27 = 42 events), plotted 
here in Galactic coordinates. Upper panel: First sector used 
for the reconstruction of the source position with the data set 
of 27 events. Red cross for the reconstructed source position; 
Middle panel: The "best sector" for the first 27 events. It 
starts from the red cross and is highlighted in orange. It con- 
tains 10 out of the 27 events. The correlation coefficient be- 
tween l/E and X' is Corr(l/E, X') ~ 0.73; Lower panel: 
The "best sector" for the first 27 events (highlighted in or- 
ange) applied to the newer 69 — 27 = 42 events. It contains 
5 events, and the correlation coefficient is ~ 0.38. Same key 
as in Fig. 0] On each panel, events located in the sectors are 
surrounded by magenta circles. See text for details on the 
"best sector". 



border point in this "best sector", the probability Pbla 
can non-negligibly vary and take for instance values as 
~ 1% or ~ 6%. This is due to the small number of 
points. The value of Pbla has to be checked in the 
future with more statistics. The order of magnitude of 
2% is compatible with the result found in the previous 



section (0.4%). 

The analysis conducted in this section is not a real 
blind analysis, since it is done a posteriori. A real blind 
analysis can start now, by fixing the best sector for the 
69 Auger events and looking at the future data. Wc 
should expect the number of events to increase in this 
sector, preferably more rapidly in average than elsewhere. 
However, the correlation coefficient Corr(X' ,1/E) will 
not necessarily increase. 



V. DISCUSSION 

The results of the two previous sections show that, in 
the Auger data, the 142 EeV event and the events in the 
Cen A region overdensity are compatible with an emis- 
sion from the Virgo cluster. 

As discussed in Ref . [37j , the Cen A region events may 
be emitted by Virgo even if they are protons. If ex- 
tragalactic magnetic fields (EGMFs) are as large as in 
Refs. [Ill UHE protons could experience deflections 
as large as several tens of degrees. 

On the contrary, if EGMFs are as low as in Refs. [13, 
Hl| , this would favour a heavy nuclei origin. Outside clus- 
ters, nuclei would not experience substantial deflections 
in the EGMF. They would be mostly deflected by the 
GMF, with such typical values. 

In the Galactic disk, the GMF is mostly parallel to the 
plane of the disk [H El- In the existing GMF models, 
the fi eld in the halo is also assumed to be parallel to the 
disk 
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In this case, nuclei from high latitude sources 
are approximately shifted along lines of equal Galactic 
longitudes. This would be consistent with the possibil- 
ity that the considered events come from Virgo. The 
142 EeV event, Cen A region and Virgo nearly have same 
longitudes. For the GMF structure, the exception is near 
the Galactic center, where a dipolar contribution may 
create a substantial component of field perpendicular to 
the Galactic plane. 

As shown in Ref. [Hf , even if there is only one source 
in the Virgo galaxy cluster, the cluster should shine as a 
whole due to substantial deflections of UHECRs in the 
magnetic fields inside. This means that the Virgo cluster 
can be considered as an extended source, basically a ~ 5° 
radius disk on the celestial sphere. Therefore, in Fig. [JJJ 
we model Virgo as a 5° radius source. 

We used the numerical code of Ref. [l4[ to deflect iron 
nuclei, with 60 EeV to 140 EeV energies, in different GMF 
models. We found that we can model the Auger data 
with several modified and reshaped versions of existing 
theoretical GMF models. For example, Figure [TU] shows 
such a modeling of the Auger data. In this figure, only 
deflections in the regular GMF are considered. If one 
adds the deflections in the turbulent field, the lower en- 
ergy events would be spread in the whole Cen A region. 
The image would then look like the considered feature 
in the Auger data. The black disk represents the Virgo 
galaxy cluster, while shaded areas show the arrival direc- 
tions of UHE iron emitted by Virgo, with a given energy. 
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FIG. 10. Image in Galactic coordinates of UHE iron nuclei 
emitted by the Virgo cluster and deflected in the regular GMF 
model which parameters are discussed in the text. Adding to 
these computations the deflections due to the turbulent GMF 
would spread the lower energy events in the whole Cen A 
region. The black disk represents the Virgo galaxy cluster. 
Shaded areas represent the arrival directions at Earth of cos- 
mic rays with given energies. Dark blue stands for 60 EeV nu- 
clei, light blue for 70 EeV, green for 80 EeV, yellow for 90 EeV, 
orange for 100 EeV, and magenta for 140 EeV. 



The dark blue region corresponds to 60 EeV iron nuclei, 
light blue to 70 EeV, green to 80 EeV, yellow to 90 EeV, 
orange to 100 EeV, and magenta to 140 EeV. 

With notations and names introduced in Ref. [l4| , the 
regular GMF model used for Figure [10] is comparable to 
a "Sun08 model" with modified parameters, and with a 
dipolar contribution similar to the "Prouza and Smida" 
model, with fiu = 30^G-kpc 3 . For the halo, we took 
B T0 = 0.8 nG, h T = 2kpc, w T , ln = 1.5 kpc, w T ,out = 
2kpc and tto = 8.5 kpc. For the disk, B = B c = 2/zG, 
r c = 5 kpc, r Q = 10 kpc, z Q = 0.2 kpc and p = 35°. Wc 
also added a —60° pitch angle in the halo. This specific 
configuration is not the only one which could lead to an 
image compatible with the data. Hence, it should not 
be considered as a prediction on the configuration of the 
GMF, or on its extension or strengths in the disk and 
the halo. It however proves that some configurations of 
the GMF are compatible with the interpretation of the 
Auger data discussed in this paper. 

The relatively low value for the correlation coefficient 
Corr(l/E,X') computed in Section [TTT1 can be notably 
due to the spread of arrival directions in the Cen A region 
due to the turbulent GMF. Such a small energy ordering 
in the image is consistent with what one can expect in 
case of heavy nuclei sources. 

If Virgo will be confirmed in the future to be a UHE 
nuclei source, it will put strong constraints on the Galac- 
tic and extragalactic magnetic fields. First, deflections in 
the EGMF would be small compared to the deflections 
in the GMF (except in the case of proton primaries, dis- 
cussed above). Second, the shift of the 142 EeV event 
would allow to give an immediate estimate of typical de- 
flections of cosmic rays in the Northern halo of the GMF. 
Third, the small scattering of arrival directions of cosmic 
rays around the thin linear structure they would have 
had in the regular field alone (see for example Fig. [TTJJ) , 



would imply that the deflections in the GMF are mostly 
dominated by the regular field contribution. The deflec- 
tions in the turbulent field would be small enough not to 
destroy the image at the lowest energies. 

The confirmation that the events in the Cen A region 
have been emitted by Virgo would put additional con- 
straints on the regular GMF: 

First, we noticed in Ref. 



141 | that large dipolar or 



toroidal contributions to the GMF can make UHE nuclei 
sources at the Galactic poles invisible. The data from 
future radio experiments will enable us to have a better 
knowledge on the strength and extensions of these com- 
ponents. If the Virgo origin of the Cen A region events is 
proved, it would independently bring strong constraints 
on the maximum contributions of the dipolar and toroidal 
components. 

Second, it would also put tight constraints on the disk 
field. The lower energy part of the image is in the Cen A 
region, which is not far from two stronger field regions: 
both the Galactic center direction and the Galactic plane. 
The computations of the images of UHE iron from Virgo, 
deflected in different GMF models, show that the shape 
of the image is very sensitive to the exact GMF config- 
uration. For heavy nuclei, the influence of the disk field 
starts to be substantial in the region b < 30° — 40° , if its 
typical height extension is non-negligible (for example for 
zq > 1 kpc). The alignment of the image along constant I 
(I ~ —60° to —30°) from Virgo to the plane would enable 
one to exclude several configurations of the field. 

Third, in our modeling, pitch angles from ~ —40° to 
~ —60° in the halo reproduced well the Auger data. This 
may suggest a non-negligible pitch angle in the halo field. 

This would also lead to a better understanding and 
tighter constraints on UHECR sources. Only a few ex- 
treme astrophysical objects can accelerate particles to 
such energies [5(| [H|. Let us note that the Auger 
UHECR flux in the Cen A region is of the same order as 
the gamma ray flux from M87, measured by HEGRA (52| 
or by H.E.S.S., MAGIC and VERITAS [Hllll. Being 
able to discriminate between the scenarii of an image cre- 
ated by one source or by several different sources in the 
Virgo cluster would require much more statistics than 
currently available. 

Reference (55| adds an important constraint on the 
Cen A region overdensity. Due to their equal rigidities, 
60 - 80 EeV iron nuclei would be deflected as 2 - 3 EeV 
protons, whatever the strengths and structures of the 
EGMF and GMF are. Hence, if UHECR sources acceler- 
ate both nuclei and protons and if this Cen A overdensity 
is made of heavy nuclei, one should expect protons at a 
Z times lower energy, exactly in the same region. The 
Auger flux is however compatible with isotropy at low 
energies, 2 - 3EeV_j55|. So if these events are nuclei, 
the results of Ref. [55[ imply either that the source spec- 
trum is harder than a 1/ E 2 spectrum or that the ratio 
of accelerated protons to nuclei in the source is not more 
than one to one. The source(s) in Virgo can have a very 
hard spectrum as, for example, in the model presented 
in Ref. 15a. A fraction of the emitted nuclei are de- 
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stroyed on their way to the observer and produce lighter 
nuclei. The events in the Cen A region would correspond 
to nuclei with a Lorentz factor of T ~ 10 9 . As shown in 
figure 1 (left panel) of Ref . [42[ , such nuclei have a mean 
free path of the order of 100 Mpc, which is much larger 
than the distance to Virgo. Besides, the propagation of 
nuclei in galaxy clusters has been studied in Ref. 
The authors find that the results mostly depend on the 
source position, as well as on the strength and profile of 
the magnetic field in the cluster. Therefore, if most of 
the emitted nuclei manage to escape the cluster, only a 
small fraction will be destroyed and enhance the light nu- 
clei flux at lower energies. Let us note that for energies 
above a few times 10 20 eV, iron and intermediate nuclei 
have a mean free path smaller than the distance to Virgo 
(see figure 1 of Ref. [HI). Therefore, the maximum ac- 
celeration energy of the source should not be too high. 
Otherwise, an additional flux of protons due to the dis- 
integration of its highest energy iron nuclei could be seen 
at lower energies in the data. Thus the confirmation of a 
nuclei source in the Virgo cluster would put interesting 
constraints on acceleration mechanisms, on the composi- 
tion of particles accelerated in the source and on physical 
conditions in the Virgo cluster. 

We show below that it will be possible to confirm or 
rule it out in the future, when more experimental data 
will be available. At the same time, we present other 
possibilities which can explain the present Auger data. 

One can expect that Auger South experiment will 
triple its statistics during its lifetime. It will confirm 
if the overdensity in the Cen A region is not a statisti- 
cal fluctuation. Auger North experiment would also be 
useful to check if there are comparable features in the 
Northern hemisphere. 

If the Cen A overdensity really exists, there are two 
cases. It is either due to protons, or to nuclei. 

If the results of HiRcs on the composition of primaries 
are correct, the events in the Cen A region are protons. 
If deflections in the EGMF arc as low as in Refs. [47l.l48j. 
most of the Cen A region events should be protons emit- 
ted by the Centaurus galaxy cluster [l|. Refs [13, [HH 
argue that this explanation would be challenged by the 
lack of events in the Virgo direction. Even if the Auger 
exposure is smaller in the direction to Virgo than in the 
direction to Centaurus, Virgo is closer to our Galaxy and 
one should statistically see at least a few protons com- 
ing from its direction. The Centaurus cluster lies behind 
Cen A. Cen A has been regarded as a potential source 
of UHE protons for a long time [59M6T| . For such a com- 
position, it may be the source of two cosmic rays in this 
direction Ja3, E3, ■ 

If the results of Auger composition studies are cor- 
rect, the events in the Cen A region are nuclei. In this 
case, there arc currently three main explanations. Ei- 
ther the UHECR deflections in the GMF and EGMF 
are large enough to prevent the identification of nuclei 
sources, or at least one source can be detected. In the 
first case, the higher flux of UHECR in this direction 
may just be due to magnetic lensing effects. Such effects 



have been studied in a particular example of lens geom- 
etry by Ref. [63j]. They have been quantified for UHE 
iron propagated in models of the regular and turbulent 
GMF by Refs. 0, HH- In the second case, there are 
two possibilities. First, these nuclei may have been emit- 
ted by Virgo and shifted in the GMF, as studied in the 
present pap er. Second, they may have been emitted by 
Cen A [40]. However, there are arguments that Cen A 
may not be powerful enou gh t o accelerate cosmic rays to 
such extreme energies jM - roq . 

There are two requirements to confirm that Virgo is 
the source. One must prove that both the 142 EeV event 
and the Cen A region events are connected to Virgo. 

Proving the link between Virgo and the highest en- 
ergy event, should be easier than for the Cen A region 
events. If other events with E > 100 EeV come in the 
region of the 142 EeV event, and are located at places 
approximately compatible with a collective emission by 
the Virgo cluster, one could prove their Virgo origin. If 
Virgo is the source, Auger may detect such an event dur- 
ing its lifetime. It would be a hint. However, the final 
confirmation should only be given by the next generation 
of UHECR experiments. 

Checking the link between the Cen A region and Virgo 
will require data from the next generation of UHECR 
experiments, such as JEM-EUSO. A confirmation of the 
Virgo origin of the highest energy event will not auto- 
matically imply that the Cen A region events have been 
emitted by Virgo. The 142 EeV event may be a nucleus 
from Virgo and be disconnected from the Cen A region 
events which may have another source. Finding some 
events in between, with more intermediate energies would 
be particularly valuable to validate that the Cen A region 
events come from Virgo. JEM-EUSO will have one or- 
der magnitude more data. It is expected to reach in few 
years of observation an exposure of 10 6 km 2 -sr-yr [67}. It 
will detect more than 1000 events at such energies. If the 
events are nuclei from Virgo, the link between the events 
with the highest energies and the Cen A lower energy 
overdensity should become clearer, suggesting a common 
origin. 

On the contrary, if the Cen A region events are protons 
from the Centaurus cluster, JEM-EUSO must be able 
to see the first characteristic individual images of UHE 
proton sources. 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper, we have proposed a new method to 
search for images of UHE heavy nuclei sources in the 
data, on top of background. 

We have pointed out that for some GMF configura- 
tions, and for some source positions on the sky, one can 
still have roughly "enlarged proton-like" images at ener- 
gies E > 60 EeV, even for iron primaries. In this case, 
one can detect the sources if they are bright enough and 
reconstruct their positions on the celestial sphere with 
the method we presented in Section [II] 
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Detecting a source in this way would however not al- 
ways be possible, because images of iron sources can often 
exhibit more complicated patterns. In general, a much 
better knowledge on the GMF (and EGMF) than cur- 
rently available is necessary to detect UHE heavy nuclei 
sources and reconstruct their positions on the sky. Future 
radio telescopes, such as SKA and its precursor LOFAR, 
will increase the number of rotation measures on the sky 
by a few orders of magnitude. This will, for example, 
enable us to know the geometry of the regular GMF in 
the halo and in the disk, as well as the turbulent GMF 
strength and properties [H, [69| . 

In Section Uni we checked if one could find such a case 
in the data recently published by the Pierre Auger Col- 
laboration [l[ . We found that the Cen A region overden- 
sity and the 142 EeV event may be the image in UHE nu- 
clei of the Virgo cluster, deflected by the GMF. With our 
method, the associated source position is reconstructed 
near the Virgo cluster, at only ~ 8.5° from M87. This 
indicates that these events are compatible with a com- 
mon origin from the Virgo cluster. The probability to 
have such a feature in some random background and re- 
construct the source at less than 10° from M87 is about 
3 x 10 -5 . If one assumes that the Cen A region overden- 
sity is due to another reason, and that the 142 EeV event 
appeared by chance near this region, the probability to 
reconstruct with our method the source at less than 10° 
from M87 is ~ 0.4%. In Section IIV1 we performed a 
"blind-like" analysis, by dividing the Auger data set in 
two parts: the first 27 Auger events and the 69 — 27 = 42 



remaining events. We determined for the 27 first events 
the "sector" on the celestial sphere for which the prob- 
ability to reproduce the data with some random back- 
ground was the lowest. We fixed it and analyzed the 
69 — 27 = 42 newer events with this sector. The proba- 
bility to have by chance the signal detected in the sector 
for the second data set is ~ 2%. 

If future data confirm that the feature discussed in this 
paper is due to UHE nuclei from Virgo, it would lead 
to significant improvements in our knowledge both on 
the cosmic magnetic fields and on the UHECR accelera- 
tion mechanisms. It would imply that deflections in the 
extragalactic magnetic fields are negligible compared to 
deflections in the Galactic magnetic field. Moreover, de- 
flections would be dominated by the regular GMF, which 
structure and strength would be better constrained. 

Thus, we have presented here both a new method to 
look for ultra-high energy nuclei sources on top of back- 
ground, and a new and consistent way to interpret the 
Auger data. We have found that one (or several) ultra- 
high energy nuclei source(s) in the Virgo galaxy cluster 
could explain both the composition and the anisotropy 
in the Auger data. However, a better knowledge of the 
Galactic magnetic field than currently available, or more 
UHECR data are still needed to confirm or rule out this 
possibility. 
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